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HIGHLIGHTS 


•  A  macro-homogeneous  model  is  developed  to  operate  Li-air  battery  in  ambient  air. 

•  The  Li-air  model  is  combined  with  an  oxygen-selective  membrane. 

•  The  model  shows  the  battery  deterioration  when  using  air. 

•  The  Li-air  performance  increases  with  integration  of  oxygen-selective  membrane. 


ARTICLE 


N  F  O 


A  B  S  T  R 


C  T 


Article  history: 

Received  17  July  2013 
Received  in  revised  form 
1  October  2013 
Accepted  24  October  2013 
Available  online  5  November  2013 


Keywords: 

Rechargeable 
Li-air  battery 

Macro-homogeneous  model 
Ambient  air  operation 
Discharge  products 
Oxygen-selective  membrane 


A  macro-homogeneous  model  has  been  developed  to  evaluate  the  impact  of  replacing  pure  oxygen  with 
ambient  air  on  the  performance  of  a  rechargeable  non-aqueous  Li-air  battery.  The  model  exhibits  a 
significant  reduction  in  discharge  capacity,  e.g.  from  1240  to  226  mAh  g^bon  at  0.05  mA  cnr2  when  using 
ambient  air  rather  than  pure  oxygen.  The  model  correlates  the  relationship  between  the  performance 
and  electrolyte  decomposition  and  formation  of  discharge  products  (such  as  Li202  and  Li2C03)  under 
ambient  air  conditions.  The  model  predicts  a  great  benefit  of  using  an  oxygen-selective  membrane  on 
increasing  capacity.  The  results  indicate  a  good  agreement  between  the  experimental  data  and  the 
model. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Several  times  higher  specific  energy  than  Li-ion  batteries  makes 
the  rechargeable  Li-air  a  candidate  of  new  generation  of  energy 
storage  devices  [1—4],  Since  the  first  report  of  a  non-aqueous 
electrolyte  Li-air  battery  in  1996  [5],  especially  after  a  break¬ 
through  in  cycle  life  by  Bruce’s  group  [6],  the  rechargeable  Li-air 
battery  has  been  investigated  intensively  [6—11],  However,  many 
science  and  technical  challenges  have  to  be  overcome  to  realise  the 
potential  of  this  cutting-edge  technology.  A  key  area  is  to  gain  in¬ 
sights  into  chemical/electrochemical  processes  that  take  place  in¬ 
side  the  Li-air  battery  via  an  effective  modelling. 
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An  aprotic  Li-air  battery  contains  a  metal  lithium  anode,  a  solid 
separator  and  a  porous  air  electrode  filled  with  a  non-aqueous  Li+ 
electrolyte.  The  fundamental  electrochemical  reactions  are  shown 
below: 

Anode  :  Li(metal)  <->Li+  +  e  (Ere{  =  0  v)  (I) 

Cathode  :  2Li+  +  02  +  2e~  <-  (Li202)solid(Erev  =  2.96  V  vs  Li) 

(II) 

Here,  Eq.  (II)  is  an  ideal  electrochemical  reaction,  desired  to 
make  a  truly  rechargeable  Li-air  battery  [6,12,13].  In  practice,  there 
are  other  side  reactions,  forming  detrimental  products,  such  as 
Li2C03  [6,12,13], 

However,  almost  all  the  current  studies  of  Li-air  batteries  are  run 
in  pure  oxygen  and  controlled  a  dried  atmosphere  in  a  glove  box  to 
minimise  the  contaminated  substances  from  ambient  air,  thus  this 
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battery  can  perform  in  high  rate  capability,  owing  to  high  oxygen 
concentration,  and  maintain  the  long-term  operation.  Our  recent 
model  studies  also  simulated  the  Li-air  battery  including  the  elec¬ 
trolyte  degradation  and  operating  in  pure  oxygen  [14,15],  To  suc¬ 
cess  in  making  a  Li-air  battery  for  use  in  practical  applications,  one 
critical  problem  to  be  solved  is  operating  Li-air  batteries  in  ambient 
air  environment  [16,17],  There  are  several  challenges  when  Li-air 
batteries  are  applied  with  air  condition  comparing  to  those  of 
pure  oxygen  feeding.  The  first  is  unavoidable  moisture  (about  1%  in 
volume)  in  surrounding  air,  which  may  penetrate  into  the  cell 
system  together  with  oxygen.  The  present  of  moisture  can  corrode 
the  metallic  lithium  anode  due  to  the  hydrolysis  reaction  with 
highly  reactive  lithium  as  shown  in  reaction  III  [17],  This  results  in 
the  fast  battery  failure  and  causes  a  serious  safety  issues. 

2Li  +  2H20<->2Li0H  +  H2  (III) 

The  second  is  the  insufficient  concentration  of  oxygen  due  to  its 
low  partial  pressure  in  atmosphere,  leading  to  limit  Li-air  batteries 
from  high  discharge  rates  because  of  small  oxygen  solubility  in  the 
electrolyte.  The  third  problem  is  that  the  small  amount  of  carbon 
dioxide  (CO2)  from  air  feeding  may  react  with  superoxide  anions, 
which  is  formed  during  the  initial  oxygen  reduction  on  discharge 
process,  to  generate  carbonate  species,  leading  to  the  deposition  of 
these  by-products  on  the  cathode  surface  [12,18,19],  It  is  worth 
noting  that  CO2  could  be  considered  as  the  active  material  to  form 
the  Li2CC>3  and  lithium  alkyl  carbonates  instead  of  desired  produce 
Li202  [12,19,20],  Furthermore,  these  side  reactions  lead  to  contin¬ 
uous  and  irreversible  consumption  of  electrolytes  and  thus  the  Li- 
air  batteries  cannot  maintain  their  sustainable  ability  during 
charge/discharge  cycles. 

To  minimise  the  battery  problem  when  using  air,  Li-air  batteries 
under  current  research  are  mainly  operated  in  pure  oxygen  [5-7,9], 
There  is  a  small  volume  of  published  works  demonstrating  the 
behaviour  of  Li-air  battery  with  air  feeding  [11,16,17,21],  These 
studies  solved  the  contaminated  gases  problems  by  using  an 
oxygen-selective  water  barrier  membrane  covering  the  outer  sur¬ 
face  of  the  cathode  to  prevent  moisture  and  permeate  oxygen  into 
the  porous  cathode  at  the  same  time.  Zhang  et  al.  developed  an 
oxygen-selective  immobilised  liquid  membrane  for  a  non-aqueous 
Li-air  battery  operated  in  ambient  air  with  20-30%  relative  hu¬ 
midity  [17],  The  membranes  were  easily  prepared  by  soaking  high 
viscosity  silicone  oil  into  porous  metal  or  Teflon  substrates.  A  Li-air 
integrated  with  these  membranes  can  be  operated  in  ambient  air 
for  16.3  days  with  a  specific  capacity  of  789  mAh  g_1  and  specific 
energy  of  2182  Wh  kg-1  based  on  the  weight  of  carbon.  The  same 
group  also  studied  the  hydrophobic  zeolite  membrane  and  poly(- 
tetrafluoroethylene)  (PTFE)  membrane  as  oxygen-selective  water 
barriers  [16],  The  latter  protected  a  Li-air  against  moisture  and 
supplied  oxygen  for  21  day  with  a  specific  capacity  of  1022  mAh  g-1 
and  specific  energy  of 2792  Wh  kg-1  based  on  the  weight  of  carbon. 

Zhang  et  al.  investigated  the  ambient  operation  of  non-aqueous 
Li-air  batteries  integrating  with  a  heat-sealable  polymer  membrane 
to  serve  as  both  an  oxygen-diffusion  membrane  and  moisture 
barrier  [11],  The  membrane  could  also  reduce  the  evaporation  of 
electrolytes  during  battery  operation.  The  Li-air  battery  with  this 
membrane  demonstrated  the  discharge  capability  in  ambient  air  for 
more  than  one  month  with  a  specific  energy  of 362  Wh  kg-1  based 
on  the  total  weight  of  the  battery  including  its  packaging.  However, 
all  the  previously  mentioned  ambient  air  works  only  studied  the  Li- 
air  performances  on  discharge  phase  without  showing  those  on 
charging  or  cycling  process,  which  could  provide  more  important 
data  for  battery  stability  than  only  single  discharge. 

In  this  paper  our  previous  Li-air  model  with  electrolyte  degra¬ 
dation  is  modified  to  operate  the  Li-air  battery  in  ambient  air 


environment  [14,15],  which  severely  damages  Li-air  performance 
and  is  still  a  critical  problem  to  be  solved  before  the  Li-air  battery 
can  be  used  for  practical  application.  The  two  species  of  oxygen  and 
CO2  are  considered  in  the  model  as  air  feeding  condition  with  the 
exception  of  moisture  or  using  dried  air  in  the  model.  Although  this 
model  does  not  include  all  the  effects  that  could  occur  in  the  Li-air 
battery  with  ambient  air  operation,  at  least  we  carefully  consider 
the  key  mechanisms  covering  the  main  impacts  of  using  air  that 
directly  affect  the  porous  cathode  behaviour  during  cell  operation, 
leading  to  detrimental  cell  performance.  Therefore,  mathematical 
model  for  Li-air  battery  with  air  feeding  condition  can  be  used  to 
identify  cell-limiting  mechanisms  and  reduce  the  time-consuming 
work  compared  to  experiment  aspect.  Moreover,  it  also  avoids  the 
serious  safety  problems  that  could  happen  when  the  ingress  of 
moisture  reacts  with  lithium  metal  anode.  This  model  can  be  used 
to  describe  the  behaviour  of  Li-air  batteries  in  ambient  air  condition 
as  well  as  to  optimise  the  performance  and  structure  of  these 
battery  electrodes. 

2.  Theoretical  mechanism  analysis 

Atypical  Li-air  battery,  shown  in  Fig.  la,  contains  a  lithium  metal 
anode,  a  separator  containing  electrolyte,  and  a  porous  carbon  or 
catalyst-loaded  carbon  air  electrode  filled  with  an  organic  elec¬ 
trolyte  comprising  a  dissolved  lithium  salt  in  an  aprotic  solvent. 
During  discharge  the  lithium  metal  anode  oxidises  to  Li+  and 
electrons  conduct  through  the  external  circuit,  while  Li+  transports 
towards  the  porous  cathode.  Oxygen  is  reduced  at  the  active  surface 
with  Li+,  thus  leading  to  the  desired  discharge  products  of  Li202  and 
the  by-product  of  L^COs  or  lithium  alkyl  carbonates  resulting  from 
the  electrolyte  decomposition  [12,13,19],  These  products  influenced 
the  Li-air  performance  and  cannot  completely  remove  during  bat¬ 
tery  cycling.  To  simplify  our  simulation,  the  model  assumed  that 
Li202  was  the  main  discharge  product  depositing  inside  the  porous 
cathode  (Eq.  (II))  and  the  irreversible  L^CCb  by-product  coexisting 
with  Li202  when  using  non-aqueous  electrolytes.  This  section  de¬ 
scribes  the  two  main  mechanisms  that  occur  inside  the  Li-air  bat¬ 
tery  and  were  applied  to  the  model. 

2.1.  Effect  of  using  ambient  air  condition 

The  critical  problems  when  using  air  as  a  feeding  reactant  for 
the  Li-air  battery  are  the  low  oxygen  solubility  in  the  electrolyte 
and  the  CO2  gas  diffusing  into  the  battery  with  oxygen.  To  inves¬ 
tigate  these  effects,  the  Li-air  model  was  changed  the  feeding 
condition  from  pure  oxygen  to  ambient  air  at  the  porous  cathode 
entrance.  Considering  the  amount  of  oxygen  and  CO2  content  in  the 
atmosphere  [22]  (78%  N2, 21%  O2, 0.035%  CO2),  the  concentration  of 
oxygen  and  CO2  in  the  electrolyte  can  be  determined  in  term  of 
Henry’s  law  which  states  that 

cg  =  ^gPg  (1) 

where  cg  is  the  concentration  of  gasses  in  the  electrolyte,  Hg  is  the 
Henry’s  law  constant  which  depends  on  the  electrolyte  and  tem¬ 
perature  used  in  the  Li-air  battery,  and  pg  is  the  partial  pressure  of 
the  gas  which  depends  on  mole  fraction  of  each  specie  in  the  at¬ 
mosphere.  Henry’s  law  is  correct  to  describe  the  solubility  of  gas  for 
low  concentrations  and  low  partial  pressures.  The  concentrations  of 
oxygen  and  CO2  calculated  from  Henry’s  law  can  be  compared  in 
Table  1.  These  concentrations  are  applied  in  the  Li-air  model  to 
represent  the  air  feeding  condition,  unless  specified  otherwise.  It 
can  be  seen  that  although  the  CO2  composition  is  not  much  in  the 
atmosphere,  its  Henry’s  law  constant  is  almost  two  times  higher 
than  that  of  oxygen,  i.e.  the  solubility  of  CO2  in  non-aqueous 
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Fig.  1.  Schematic  computation  domain  of  a  Li-air  battery  during  discharge  operation,  (a)  Li-air  battery  operated  with  ambient  air  feeding  and  the  inset  demonstrates  the  discharge 
products  formation  of  Li202  and  Li2C03  covering  on  the  porous  carbon  surface,  (b)  Li-air  battery  protected  by  an  oxygen-selective  membrane  at  the  cathode. 


solvents  is  higher  than  oxygen  [23].  This  high  solubility  of  CO2  into 
the  electrolyte  could  affect  the  Li-air  performance  due  to  the  for¬ 
mation  of  Li2CC>3. 

Hence,  unlike  our  previous  model  using  only  pure  oxygen  as 
feeding  condition,  the  present  model  applied  the  air  condition  with 
low  solubility  of  oxygen  and  CO2.  These  concentrations  in  the 
porous  cathode  can  affect  the  Li-air  performance  during  battery 
cycling.  Nitrogen  also  accesses  through  the  cathode  but  has  a  little 
or  no  effect  with  lithium-based  electrolyte  on  the  Li-air  battery 
performance  [24], 

2.2.  Electrolyte  degradation 

In  order  to  succeed  in  the  application  of  non-aqueous  Li-air 
batteries,  finding  of  electrolytes  with  high  stability  during  cell 
operation  is  a  prerequisite  for  long-cycled  life  of  Li-air  battery, 
especially  in  an  oxygen-rich  environment  and  catalytic  condition 
[25],  Since  the  Li-air  battery  have  been  developed,  there  is  still 
no  proper  electrolyte  that  can  meet  all  of  these  requirements:  1) 


high  stability  during  cell  operation  [12,18,26,27],  2)  high  boiling 
point  (low  evaporation),  and  3)  high  oxygen  solubility  and 
diffusivity  as  well  as  less  sensitive  to  moisture  [28],  At  present, 
various  non-aqueous  electrolytes  have  been  widely  examined 
and  applied  in  Li-air  batteries  based  on  different  types  of  sol¬ 
vents,  e.g.  carbonates,  ethers,  sulfoxides,  nitriles  and  ionic  liquid 
[29],  However,  Li2C03  formation  from  decomposition  of  elec¬ 
trolyte  by  Li— O2  intermediates  and  products  was  also  clearly 
detected  inside  the  electrode  during  discharge  from  all  electro¬ 
lytes  dissolved  in  different  solvents  [29],  Hence,  the  electrolyte 
degradation  mechanisms  proposed  in  previous  work  [15]  are 
also  included  in  the  present  model  and  reinstated  as  described 
below. 

2.2.1.  U2CO3  formation 

202  +  2e-^202£°  =  3.18  (la) 

202  +  2C02^C20g_  +  02  RDS  (lb) 
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Table  1 

Henry’s  constants  and  calculated  concentrations  in  non-aqueous  at  25  °C. 

Species  Henry’s  law  constant  Concentration  Ref. 

(mol  m~3  atm-1)  (molm-3) 

02  2.95  0.6182  Calculated  from  Ref.  [47] 

C02  4.59  0.0043  [48] 


C20 +  4Li+  +  202  — >2Li2C03  +  202  (lc) 

Overall  4Li+  +  2C02+202  +2e  ^2Li2C03  +  02  £°  =  7.19 

(IV) 

2.2.2.  Solvent  degradation 


02  +  e-^0;T  E°  =  3.18  (2a) 

02  +  [Solvent]  ->  [Solvent]  (2b) 

[Solvent]-  +  02^C02  +  H20  +  HC02Li  +  CH3C02Li  (2c) 

Overall  202  +  [Solvent]  +  e  ->■ C02  +  H20  +  HC02Li  +  CH3C02Li 

(V) 

It  is  worth  noting  that  the  exact  details  of  the  reaction  routes  for 
the  electrolyte  degradation  to  form  Li2C03  can  be  complicated 
regarding  several  intermediates  during  cell  operation.  Then,  the 
proposed  mechanisms  above  are  a  possibility  to  describe  the  for¬ 
mation  reaction  of  Li2C03  as  reported  previously  [19,30,31], 

With  the  use  of  air-feeding  condition  instead  of  favourite  pure 
oxygen  to  the  porous  carbon  cathode,  the  C02  in  this  model  cri¬ 
terion  could  be  generated  from  2  routes,  one  from  solvent  or 
electrolyte  degradation  (Eq.  (V))  and  the  other  from  the  atmo¬ 
spheric  air.  Both  Li202  and  Li2C03  can  be  produced  as  discharged 
products  and  are  usually  insoluble  in  the  cell  electrolyte.  As  a 
result,  the  repeated  depositing  film  of  different  lithium  salts  over 
the  carbon  surface  after  each  discharge  and  charge  step  affects  the 
reaction  mechanisms  at  the  active  area  and  electrolyte  interfaces, 
and  so  the  species  transport  inside  the  pores  leading  to  cell 
voltage  loss  [32], 

3.  Battery  model 

3.1.  Model  description  and  assumptions 

This  work  investigates  the  impacts  of  using  air  on  a  one¬ 
dimensional  model  for  a  prismatic  single  cell  of  a  Li-air  battery 
which  constitutes  of  a  thin  lithium  metal  as  negative  electrode,  an 
anode  protective  layer  (APL),  a  porous  separator  soaked  with  non- 
aqueous  electrolyte,  and  a  porous  carbon  air  cathode  served  as 
reaction  site  and  accommodation  for  discharge  products  as  shown 
in  Fig.  la.  Current  collectors  are  placed  at  the  back  of  each  elec¬ 
trode.  The  electrolytic  solution  is  considered  as  a  concentrated 
binary  electrolyte  to  describe  the  motion  for  of  each  species  in  the 
systems.  Due  to  the  larger  thickness  of  porous  cathode  than  both 
of  APL  and  separator,  this  area  represents  as  the  main  contribution 
to  influence  the  cell  performance,  thus  we  are  particularly  inter¬ 
ested  in  analysing  the  effect  of  different  parameters  on  the  porous 
cathode  for  the  Li-air  simulation.  The  formation  of  solid  products 
inside  the  porous  cathode  is  based  on  a  macro-homogeneous 
porous  model  defining  the  electrode  by  its  porosity  which  is 


initially  uniform  but  which  dynamically  changes  as  discharge 
occurs. 

As  the  Li-air  battery  system  is  complex  with  the  various  re¬ 
actions  and  mass  transport  species  along  the  entire  cell,  the  build¬ 
up  of  Li202  and  Li2C03,  and  the  change  of  porosity  and  interfacial 
surface  area  vary  with  time  and  space,  several  model  assumptions 
were  adopted  to  support  the  calculation  as  follow: 

•  The  Li-air  battery  is  operated  in  isothermal  conditions  so  that 
the  thermal  effects  are  not  considered 

•  The  lithium  salts  of  Li202  and  Li2C03  are  the  main  discharge 
products  which  only  occur  and  deposit  inside  the  porous  cathode. 

•  The  electrolytes  used  in  Li-air  batteries  are  assumed  a  binary 
monovalent  electrolyte  which  consists  of  a  single  salt  in  a  ho¬ 
mogeneous  organic  solvent  mixture. 

•  The  electrolyte  behaviour  is  based  on  concentrated  solution 
theory  to  simulate  the  Li+  diffusion. 

•  The  pores  in  the  cathode  are  filled  of  liquid  phase  electrolyte 
such  as  a  solution  of  lithium  hexafluorophosphate  (LiPFe)  or 
lithium  bistrifluoromethanesulfonamide  (LiTFSI)  in  a  non- 
aqueous  solvent. 

•  The  dried  air  feeding  is  assumed  to  dissolve  as  the  concentration 
in  the  organic  electrolyte  with  a  saturated  initial  concentration 
so  that  no  gas  phase  is  occurred  in  the  system. 

•  Convection  for  mass  transport  is  negligible  inside  the  cell. 

3.2.  Governing  equations 

The  Li-air  model  applied  in  this  work  is  based  on  previously 
developed  literature  in  micro-macro  homogeneous  model  for 
electrochemical  cell  system  and  Li-air  battery  [14,15],  This  section 
provides  the  governing  equations  which  describe  conservation  of 
mass  and  current,  species  transport,  and  reaction  kinetic  in  the 
porous  cathode  and  separator  to  clarify  the  mechanism  inside  the 
Li-air  battery  with  ambient  air  operation. 


3.2.1.  Transport  of  species 

As  the  model  is  based  on  the  macroscopic  theory  of  porous 
electrode  which  considers  the  solution  and  solid  matrix  phases  as 
superimposed  continuum  [33,34],  Based  on  this  approach,  a  ma¬ 
terial  balance  equation  for  species  i  transport  in  the  Li-air  electro¬ 
lyte  can  be  expressed  as 


a(£Ci) 

at 


(2) 


where  c,-  is  the  bulk  concentration  of  species  i  in  the  solution  phase 
which  is  averaged  over  the  volume  of  the  solution  in  the  pores,  e  is 
the  porosity  of  the  electrode  which  is  the  electrolyte  space  in  the 
matrix  phase,  JV,  is  the  molar  flux  of  species  i  in  the  porous  solution 
averaged  over  the  cross  sectional  area  of  the  electrode,  and  r;  is  the 
volumetric  production  rate  of  species  i  from  the  solid  phase  (elec¬ 
trode  material)  to  solution  phase  (electrolyte  in  the  porous)  within 
the  porous  electrode. 

The  concentration  of  lithium  salt  electrolyte  is  the  same  as 
concentration  of  Li+  due  to  the  binary  electrolyte  assumption. 
Without  convection,  the  diffusion  and  migration  fluxes  equations 
for  mass  transfer  of  Li+  and  all  of  the  species  transporting  in  the 
porous  cathode  can  be  expressed  as 


JVLi  =  -DLi,effVCLi+Ip-  (3) 

Ni  =  — DjeffVCj  (4) 

where  Dyieff  and  Djeff  are  the  effective  diffusion  coefficient  of  Li+ 
and  species  i,  respectively,  t+  is  the  transference  number  of  Li+,  F  is 
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Faraday’s  constant  which  is  equal  to  96,485  C  mol-1,  and  i-2  is  the 
current  density  in  the  solution  phase  or  electrolyte  current  density 
which  can  be  defined  by  the  gradient  of  the  potential  in  a  1 : 1  binary 
concentrated  electrolyte  solution  as  [33,35] 


h 


-KeffVfo  - 


2RTKef[ 

F 


(4  - 1) 


ain f) 

dlncu) 


Vlncu 


(5) 


where  Keff  is  the  effective  conductivity  of  the  electrolyte, <j>2  is  the 
electrolyte  potential  (electric  potential  of  Li+),  R  is  the  universal  gas 
constant  which  is  equal  to  8.3143  J  mol-1  K-1,  T  is  the  cell  tem¬ 
perature  in  Kelvin,  and  /  is  the  activity  coefficient  of  LiPF6  salt. 

In  the  solid  matrix  phase,  the  movement  of  electron  is  governed 
by  Ohm’s  law  which  evaluates  the  electric  potential  variation  or 
potential  of  electron,  as  follow 


coefficient  of  the  species  i,  n  is  the  number  of  electron  transferred  in 
the  reaction,  a  is  the  specific  interfacial  area  of  the  pore  per  unit 
volume  of  the  total  electrode,  and  jm  is  local  transfer  current  density 
between  electrode  and  electrolyte  interface  of  each  reaction  at  the 
cathode.  The  value  of  s,-,  z,-,  and  n  can  be  defined  by  matching  with 
an  individual  electrode  reaction  using  the  general  form  of  Eq.  (12), 
for  example,  the  value  of  su,  Zu,  and  n  of  Li+  from  Eq.  (II)  are  -2, 1, 
and  2,  respectively. 

In  practical  Li-air  cell,  the  precise  reaction  routes  can  be 
complicated  regarding  several  intermediates  as  proposed  in  the 
previous  reports  [12,19],  Hence,  both  electrochemical  reactions  of 
Li202  and  Li2C03  formation  inside  porous  electrode  are  considered 
in  the  present  work.  The  superficial  production  rate  of  each  species 
(referred  to  Eq.  (2))  from  solid  phase  to  pore  solution  in  the  indi¬ 
vidual  reactions  m  is  given  by  Faraday’s  law 


i\  =  -<teffV(4  (6) 

where  creff  is  the  effective  conductivity  of  the  electron  in  the  elec¬ 
trode.  This  parameter  is  affected  by  the  volume  fraction  of  solid 
electrode  inside  the  porous  cathode. 

Due  to  the  tortuosity  of  porous  cathode,  the  effective  parameters 
of  Dj,eff,  Keff  and  (Teff  in  the  above  equations  are  affected  by  the 
change  of  solid  volume  fraction  and  porosity  (assigned  as  variable  e 
in  the  model)  inside  the  porous  medium.  These  parameters  are 
applied  for  only  the  porous  cathode  region  (others  are  not  porosity) 
and  are  corrected  to  account  for  the  porosity  effect  using  Brugge- 
man  correlation  [36] 

Du, eft  =  c15DLi  (7) 

Dieff  =  e1'5Di  (8) 

Keff  =  e1'5*  O) 

<teff  =  (l-e)1'5*  (10) 


3.3.  Rate  expression  at  cathode 

The  actual  reaction  paths  and  mechanisms  for  the  discharge 
products  are  not  available  and  quite  complex  involving  various 
intermediates.  Hence,  to  describe  the  electrochemical  kinetic 
expression  for  the  porous  cathode,  the  model  adopts  the  kinetic 
expression  based  on  Eq.  (II)  for  Li2C>2  formation  and  on  Eq.  (V)  for 
L^CCh  formation. 

3.3.1.  U2O2  formation 

For  electrochemical  reaction  of  Li2C>2  at  the  cathode,  a  modified 
version  of  the  Butler— Volmer  equation  is  applied  in  the  model 
using  two  rate  coefficients.  The  reaction  for  Li2C>2  formation  pre¬ 
sented  in  Eq.  (II)  also  depends  on  the  concentration  of  Li+  and 
oxygen  for  discharge  and  the  concentration  of  Li202  during  charge 
as  in  the  following  equation 


where  Du,  D,-,  k  and  a  are  the  diffusion  coefficient  of  the  Li+  and  each 
specie  in  electrolyte  and  the  conductivity  of  electrolyte  and  elec¬ 
tron  in  the  cathode,  respectively. 

3.2.2.  Conservation  of  charge 

For  the  porous  electrode  theory,  the  charge  conservation  for  the 
matrix  and  solution  phases  would  require  the  divergence  of  the 
total  current  density  to  be  zero  defined  by 

Vi!  +  Vi2  =  0  (11) 

During  discharge  or  charge,  the  electrochemical  reactions 

occurring  at  the  electrode/electrolyte  interface  (charges  transfer 

reaction)  are  expressed  for  individual  reactions  according  to  the 
conventionally  general  formula  of  the  form 

Mf^ne-  (12) 

The  charge  transfer  from  solid  phase  to  electrolyte  phase  per 
unit  volume  of  electrode  (V  i2)  is  related  to  the  individual  average 
transfer  current  density  occurred  at  the  cathode  given  by 

Vi2  =  £q/m  (13) 

This  equation  states  that  the  transfer  current  per  unit  electrode 
volume  is  equivalent  to  the  electrode  chemical  reaction  rate  where 
the  Mi  is  a  species  symbol  participating  in  the  electrochemical  re¬ 
action,  z;  and  Si  are  the  charge  number  and  the  stoichiometric 


k,(cU202)exp [(1  Rj)n%m]  - 4(cLii)2(c0Jexp 

(15) 

Vm  =  <h  ~  02  -  -  £m  (16) 

A0film  =  lc4lmfc’s  (17) 

where  k3  and  kc  are  the  anodic  and  cathodic  rate  constant, 
respectively,  /?  is  the  symmetry  factor  equal  to  0.5,  ijm  is  surface  or 
activated  overpotential  for  individual  reaction,  m,  at  the  cathode, 
and  Rf,|m  are  the  voltage  drop  and  the  electrical  resistivity 
across  U2O2  film  formation,  respectively,  es  is  the  volume  fraction  of 
solid  formation  of  discharge  products  of  IA2O2  and  Li2C03,  and  E°,  is 
the  theoretical  open-circuit  potential  for  each  reaction. 

3.3.2.  U2CO3  formation 

From  our  previous  work  [15],  we  proposed  2  mechanisms  for 
Li2CC>3  formation,  one  generated  from  reaction  between  superoxide 
and  CO2  (Eq.  (IV))  and  the  other  from  the  C02  reduction  reaction. 
Because  preliminary  simulated  results  showed  that  the  Li2C03 
formation  from  the  latter  demonstrated  the  insignificant  value 
compared  to  the  former,  this  model  includes  only  one  Li2C03  for¬ 
mation  mechanism.  As  explained  before,  the  Li2C03  formation  is 
one  of  the  discharge  by-products  coexisting  with  I22O2.  Thus,  the 
decomposition  of  electrolyte,  which  initially  forms  C02  and  finally 


U.  Sahapatsombut  et  at  /  Journal  of  Power  Sources  249  (2014)  418-430 


423 


Table  2 

Parameters  used  in  simulation. 


Parameter 

Cell  properties 

Thickness  of  APL 

Thickness  of  separator 

Thickness  of  porous  positive  electrode 

Conductivity  of  positive  electrode 

Porosity 

Specific  interfacial  area  of  cathode 

Electrical  resistivity  across  Li202  film  formation 

Electrolyte  properties 

Electrolyte  concentration 

Li+  diffusion  coefficient3 

Oxygen  diffusion  coefficient 

Superoxide  diffusion  coefficient 

Carbon  dioxide  diffusion  coefficient 

Conductivity  of  Li+  in  electrolyte 

Transference  number  of  Li+a 

dlnfltflncu? 

Kinetic  parameters 

Reaction  rate  coefficient  anodic  current 
Reaction  rate  coefficient  cathodic  current 
Reaction  rate  coefficient  for  OJ  formation 
Reaction  rate  coefficient  for  Li2C03  formation 
Reaction  rate  constant  of  C02  formation 
Exchange  current  density  for  anode 
Symmetry  factor 
General  parameter 

Mass  density  of  Lithium  peroxide  (Li202) 

Mass  density  of  electrolyte  solution  (LiPF6) 
Mass  density  of  Lithium  carbonate  (Li2C03) 
Mass  density  of  carbon 
Operating  temperature 


5  x  1(T8 
5  x  10-5 
7.5  x  itr4 
10 
0.73 

3.75  x  106 
50 

1000 

2.11  x  10~9 
7  x  1010 
9  x  1010 
1  x  10~9 
1.085 
0.2594 
—1.03 


Unit 


kg  m  3 
kg  m~3 
kgm-3 
kgm-3 


K 


La 

Lc 


Rfilm 


[11] 

[11] 

[11] 

[49] 

[50] 

Calculated 

[51] 


ka  Assumed 

kc  Assumed 

kia  Assumed 

kit 

kv  [56] 


(3 


[49] 


/>u2o2 

PUPF, 

/>U2C03 

Pc  [57] 


generates  Li2C03  as  described  above,  occurs  during  Li-air  operation. 
Some  of  the  elementary  steps  in  the  electrolyte  degradation  se¬ 
quences  that  lead  to  these  by-products  may  be  irreversible  and 
non-electrochemical,  so  that  the  overall  kinetic  expression  is  very 
complex.  Therefore,  we  use  the  kinetics  for  Li2C03  formation  on  Eq. 
(IV)  above  based  on  the  published  kinetic  data. 

First  the  superoxide  radical  anion  (O2 )  that  is  initially  formed 
(Eq.  (la))  during  Li-air  discharge  as  evidenced  in  previous  study 
[37],  attacks  CO2  which  is  generated  from  solvent  decomposition 
(from  Eq.  (V))  to  finally  form  I22CO3  with  the  presence  of  Li+  as 
follow  reaction 

in  =  FMcoJ  [  -  exp(  ^"F??m)j  (18) 

rlb  =  klb(c02)  (CC02)  (19) 

where  kla  and  ktb  are  the  rate  constant  for  the  electrochemical 
reaction  to  form  O2  and  chemical  reaction  to  generate  L^COa, 
respectively.  We  use  the  Tafel  form  in  Eq.  (18)  rather  than  the 
Butler— Volmer  because  the  large  kinetic  overpotential  during  cell 
discharge  puts  the  reaction  in  the  Tafel  region  and  considers  only 
discharge  (irreversible  for  Of  formation).  It  has  been  demonstrated 
that  the  chemical  reaction  in  Eq.  (19)  is  found  to  be  first-order  with 
to  respect  to  both  O2  with  CO2  reactants  and  is  the  rate  deter¬ 
mining  step  (RDS)  [31],  Hence,  the  other  reaction  is  considered  as 
equilibrium  and  the  formation  of  Li2CC>3  can  be  predicted  by  using 
Eqs.  (18)  and  (19)  together. 

3.3.3.  Solvent  degradation 

The  CO2  generation  can  be  created  from  the  electrolyte  degra¬ 
dation  which  is  first  attacked  by  the  superoxide  formation  as 


described  in  Eq.  (V).  Apart  from  this  decomposition,  some  reports 
demonstrate  that  side  reactions  to  form  CO2  were  observed  at  the 
cathode  and  were  attributed  to  carbon  decomposition  during 
charge  process  [29,38],  However,  only  the  CO2  generated  from 
electrolyte  degradation  will  be  considered  here.  Addressing  this 
kinetic  expression,  because  a  detailed  mechanism  is  not  available, 
the  CO2  formation  based  on  the  overall  reaction  in  Eq.  (V),  which 
consider  the  solvent  concentration  as  constant  is: 

h  =  Fkv (c02)[-  exp ym ) ]  (20) 

where  kv  is  the  rate  constant  for  the  electrochemical  reaction  to 
form  CO2  and  the  others  are  the  same  as  described  above. 

3.4.  Rate  expression  at  anode 

The  electrochemical  reaction  rate  for  the  anode  includes  the 
oxidation  of  lithium  metal  to  soluble  Li+.  It  is  described  by  a  general 
Butler— Volmer  equation  as  follow 

h  =  *0  [exp((1  rt  nF??a) _  exp(zwrVi)}  (21) 

where  to  is  exchange  current  density  for  anode,  is  surface  or 
activated  overpotential  for  reaction  at  anode,  and  the  other 
parameter  are  as  described  above. 

3.5.  Specific  surface  area  and  porosity  change 

The  specific  area  (a)  of  the  electrode/electrolyte  interface  in  Eq. 
(22)  is  decreased  by  the  morphology  and  dynamic  change  of  the 
porosity  due  to  the  solid  discharge  products  of  IJ2O2  and  Li2CC>3. 
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These  products  are  insoluble  in  several  non-aqueous  electrolytes 
and  cover  the  active  surface  area  during  battery  discharging.  The 
variation  of  effective  local  surface  area  per  unit  volume  of  electrode 
can  be  commonly  written  by  a  geometric  relation  [39,40] 


where  es  and  e°  are  the  volume  fraction  of  discharge  products  of 
Li2C>2  and  Li2CC>3  solid,  and  initial  electrode  porosity,  respectively. 
This  empirical  equation  is  used  to  describe  the  change  in  the 
interfacial  area  for  electrochemical  reactions  that  occur  during 
discharge  because  of  the  fast  passivation  of  Li202  and  Li2C03 
covering  a  portion  of  the  active  sites  for  electrochemical  reaction 
over  the  carbon  surface  [41  ].  The  magnitude  of  exponent  p  is  a 
geometrical  factor  indicating  the  morphology  shape  of  the  solid 
products  that  cover  the  active  area.  Small  values  of  p  indicate  that 
the  flat,  plate-like  precipitate  shape  of  solid,  conversely,  large 
values  of  p  reflects  the  needle-like  solid  which  cover  small  active 
area.  In  this  model,  the  value  of  0.4  has  been  used. 

The  porosity  volume  change  of  the  carbon  electrode  will  be 
decreased  due  to  the  formation  of  insoluble  solid  products  covering 
the  catalyst  and  active  particles  as  described  in  Eq  (23).  Thus,  the 
effective  conductivities  and  diffusivities  for  both  Li+  and  all  the 
species  inside  the  porous  cathode  are  affected  by  the  porosity 
change.  This  can  be  described  by  the  Bruggeman  relationship 
(referred  to  Eqs.  (7)— (10)).  Because  the  solid  distributions  in  the 
model  are  U2O2  and  U2CO3,  we  use  the  properties  of  these  solids  for 
all  of  the  discharge  products  formed  on  the  individual  reaction  m 


Ot  V  >  .  IVlm  /m\ 

at  ~  (23) 

solid  phase  m  rm 

The  volume  fraction  of  the  discharge  solid  formation  can  be 
determined  from  the  cathode  volume  balance  as 


(24) 


where  Mm  and  pm  are  the  molecular  weight  and  the  mass  density  of 
solid  discharge  products,  respectively. 


current  density  in  the  electrolyte  phase  equals  to  zero,  and  the  flux 
of  each  species  is  zero.  At  the  cathode  electrode/separator  interface 
(x  =  Lc)  the  continuous  boundary  conditions  are  specified  for  the 
fluxes  of  all  species.  The  current  density  in  the  solid  phase  in  this 
interface  becomes  zero,  and  the  current  density  in  the  electrolyte 
then  equals  to  the  applied  discharge  current  density.  The  voltage  of 
the  cell  is  calculated  by  the  difference  between  the  electrode  po¬ 
tential  at  cathode  current  collector  and  the  electrolyte  potential  at 
the  anode  side,  Vceii  =  0i(x  =  L)  -  (j> 2(x  =  0). 

The  conservation  equations  and  the  boundary  conditions 
described  above  were  discretized  using  a  finite  element  method 
and  solved  in  one-dimensional  battery  system  by  commercial 
software  package  COMSOL  multiphysics  version  4.3.  The  COMSOL 
software  is  designed  to  solve  a  set  of  coupled  differential  and 
algebraic  equations.  The  battery  simulation  model  is  performed 
on  a  32  bit  Windows  platform  with  4GB  RAM,  and  Intel  Core  2 
Duo  2.93  GHz  processor.  The  different  transport  equations  and 
the  electrochemical  reactions  were  solved  as  time  dependent 
until  the  cell  voltage  reached  the  stop  condition.  The  solution 
was  considered  as  converged  solution  when  the  difference  be¬ 
tween  two  results  was  less  than  10-4  (relative  tolerance)  for  all 
variables. 

4.  Results  and  discussion 

The  one-dimensional  Li-air  battery  implemented  with  the 
model  equations  as  presented  in  the  previous  section  have  been 
simulated  and  solved  to  analyse  the  impact  of  using  ambient  air 
condition  which  is  a  severe  effects  on  Li-air  batteries  comparing  to 
the  case  of  pure  oxygen.  As  mentioned  before,  this  model  was  based 
on  previous  work  created  by  the  author  and  could  be  validated  by 
using  the  experimental  cycle  performance  obtained  in  our  labs,  due 
to  sufficient  rechargeable  data.  In  this  section,  first,  the  Li-air  bat¬ 
tery  is  considered  the  effects  of  using  dried  air  only  without  those 
from  electrolyte  degradation  to  distinguish  between  these  two 
mechanisms.  Then,  the  model  is  combined  the  two  effects  together 
to  investigate  the  cell  performance  in  term  of  discharge  specific 
capacity  and  IJ2CO3  accumulation  on  cycling  resulting  in  electrode 
passivation  and  capacity  fading. 


3.6.  Initial  conditions 

To  solve  the  governing  equation  for  the  battery  cycling  pro¬ 
cess,  initial  conditions  are  specified  for  all  the  species  concen¬ 
tration  inside  the  electrochemical  cell,  the  porosity,  the  specific 
interfacial  area,  and  the  cell  thickness.  These  initial  values 
applied  in  the  Li-air  battery  model  are  adopted  from  literature 
and  summarised  in  Table  2.  Before  starting  the  discharge  battery, 
the  concentration  for  each  species  is  assumed  to  uniformly 
distribute  at  all  locations  inside  the  cell  system  and  equals  to 
their  initial  concentration.  Because  the  thickness  of  the  porous 
cathode  is  much  larger  than  both  of  APL  and  separator,  then  this 
area  is  regarded  as  the  critical  region  of  the  computational 
domain  affecting  the  cell  performance. 

3.7.  Boundary  conditions 

From  Fig.  la,  a  schematic  view  of  the  model  cell  consists  of  four 
boundaries  and  three  domain  regions.  Contrary  to  previous  model, 
the  constant  for  oxygen  concentration  feeding  at  the  right  side  of 
the  cathode  (x  =  L)  can  be  estimated  form  partial  pressure  using 
Henry’s  law  constant  as  described  above.  At  the  current  collector  or 
the  back  side  of  the  cathode  electrode  (x  =  L),  the  current  density  in 
the  solid  phase  is  equal  to  the  applied  discharge  current  density,  the 


-  -  ■  •  0.1  mA  cm'2  in  Air 


Specific  capacity  (mAh  gc'arbon) 

Fig.  2.  Comparison  of  the  voltage-capacity  curve  for  a  non-aqueous  Li-air  battery 
model  in  different  feeding  conditions  between  pure  oxygen  and  ambient  air  at  two 
operating  rates  of  0.05  and  0.1  mA  cnr2.  The  electrolyte  contains  1  M  LiPF6  dissolved 
in  acetonitrile.  The  cathode  electrode  thickness  is  750  pm  with  porosity  of  0.73.  The 
cell  cycle  is  simulated  between  2.4  and  4.2  V  versus  Li/Li+  in  at  operating  temperature 
298.15  K. 
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4.1.  Li-air  performance  with  air  feeding 

To  compare  the  effect  of  using  air  into  the  Li-air  battery,  all 
parameters  applied  in  the  model  are  summarised  in  Table  2  and  the 
oxygen  and  CO2  concentrations  in  dried  air  are  applied  by  following 
the  atmospheric  condition  with  low  partial  pressure  as  presented 
in  Table  1.  This  model  section  excludes  the  effect  from  electrolyte 
degradation  and  considers  the  effect  from  air  feeding  only.  The 
voltage-capacity  results  on  1st  cycle  obtained  from  the  cell  at 
different  discharge  rates  are  presented  in  Fig.  2.  It  is  apparent  that 
at  the  high  discharge  rate  of  0.1  mA  cm-2  there  is  a  significant 
decrease  in  cell  discharge  capacity  provided  only  117  mAh  gcalbon 
when  using  air  as  the  active  reactant  unlike  the  cell  performance  in 
pure  oxygen  providing  the  first  discharge  capacity  around 
700  mAh  gcarbon  at  the  same  discharge  rate.  A  large  deterioration  of 
cell  capacity  can  be  attributed  to  the  very  low  oxygen  solubility  in 
electrolyte  at  partial  pressure  0.21  atm  (0.618  mol  m  3  compared  to 
pure  oxygen  of  3.264  mol  m  3  at  1  atm)  and  high  discharge  rate 
resulted  in  limited  oxygen  diffusion  in  the  porous  cathode.  The  CO2 
concentration  which  is  even  low  solubility  than  that  of  oxygen 
could  not  much  affect  the  cell  capacity  from  Li2C03  formation  on 
1st  cycle. 

Because  the  cell  performance  provided  the  low  capacity  at  high 
discharge  rate  and  could  not  determine  the  cell  cycling  behaviour, 
the  discharge  rate  in  the  model  was  changed  to  lower  value  as 
0.05  mA  cm-2  to  investigate  the  air  effect  compared  to  the  pure 
oxygen  at  the  same  cycling  rate  as  also  presented  in  Fig.  2.  It  is 
apparent  from  the  graph  that  although  the  battery  capacity  in¬ 
creases  on  changing  the  discharge  rate,  there  has  been  a  clear 
reduction  of  capacity  from  1240  to  226  mAh  gcarbon  when  using 
ambient  air  as  the  feeding.  What  is  also  interesting  in  this  data  is 
that  the  inlet  air  affected  not  only  the  cell  capacity  but  also  the 
discharge  voltage  decreasing  from  2.75  V  in  pure  oxygen  to  2.55  V 
in  air  at  the  same  discharge  rate.  The  cell  performance  is  depleted 
because  the  oxygen  substantially  reduces  in  diffusion  and  solubility 
farther  inside  the  porous  structure  when  a  Li-air  battery  operated 
in  ambient  condition.  Moreover,  the  cycling  profile  of  Li-air  battery 
in  Fig.  3  indicates  that  the  discharge  voltage  slightly  drops  due  to 
gradual  Li2C03  deposition  when  repeatedly  cycles  the  Li-air  battery. 
Because  the  electrolyte  degradation  effect  is  not  include  in  this 
model  section  to  clarify  the  ambient  air  effect  only,  then  the  battery 


Specific  capacity  (mAh  g^>rbon) 

Fig.  3.  Variation  of  voltage-capacity  curve  for  a  non-aqueous  Li-air  battery  using 
ambient  air  in  8  cycles  on  discharge  and  then  charge  between  2.2  and  4.2  V  versus  Li/ 
Li+  at  a  rate  of  0.05  mA  cnr2.  The  other  parameters  used  in  the  model  are  the  same  as 
described  in  Fig.  2. 


Cycle  number 


jn  of  8-cycle  rechargeable  Li- 


Fig.  4.  The  cycle  performance  in  term  of  capacity  re 

air  battery  in  ambient  air  feeding  compared  between  the  cells  with  stabilised  elec¬ 
trolyte  and  with  electrolyte  degradation  effect.  Battery  was  cycled  at  a  rate 
0.05  mA  cm-2.  The  other  parameters  used  in  the  model  are  the  same  as  described  in 
Fig.  2. 


profile  between  2nd— 8th  cycle  in  Fig.  3  shows  slightly  decrease  in 
specific  capacity  due  to  small  Li2C03  generated  from  low  CO2 
concentration  in  air.  The  combining  effects  of  air  and  electrolyte 
degradation  will  present  in  the  next  section. 

The  capacity  retention  of  Li-air  at  a  rate  of  0.05  mA  cm-2  is 
presented  in  Fig.  4,  when  the  cell  is  operated  by  using  an  air  inlet. 
The  graph  shows  that  the  capacity  retention  gradually  falls  from 
100%  in  the  1  st  cycle  (cell  discharge  capacity  of  226  mAh  gcalbon)  to 
65%  in  2nd  cycle  and  slightly  decreases  afterwards  to  58%  in  8th 
cycle.  The  significant  difference  of  capacity  retention  between  1st 
and  2nd  cycle  is  due  to  incompletely  reversible  Li2C>2  during  the 
recharging  cell  after  1st  discharge.  Thereafter,  the  cell  is  affected  by 
the  Li2C03  formation  (generated  from  Eq.  (IV))  and  demonstrates  a 
small  decrease  of  capacity  retention  during  the  cycling  due  to  little 
contaminated  CO2.  As  can  be  seen  from  Table  3,  the  volume  fraction 
of  Li2C03  collected  at  cathode/current  collector  interface  at  the  end 
of  each  discharge  cycle  slightly  increases  and  gradually  deposits 
inside  the  cathode  with  the  repeated  cycle  of  the  Li-air  battery.  This 
volume  fraction  is  also  plotted  in  Fig.  5  but  the  value  is  unnoticeable 
due  to  the  scale  when  compared  to  Li2C03  formation  from  the  Li-air 
model  including  the  electrolyte  degradation  effect,  which  will  be 
described  in  the  next  section.  Because  the  Li2C03  formation  tends 


Table  3 

Specific  discharge  capacity  and  U2CO3  formation  for  a  non-aqueous  Li-air  battery 
operated  in  ambient  air  in  8  cycles  at  a  rate  of  0.05  mA  cm-2.  The  data  were 
compared  between  the  cell  with  stabilised  electrolyte  and  with  electrolyte  degra¬ 
dation  effect. 


Cycle  Discharge  capacity  Capacity  U2CO3  volume 

(mAh  gcarbon)  retention  (%)  fraction 


Feeding  air  WithEDa  Feeding  ail 
226.65  250.47  100.00 

146.84  157.70  64.79 

143.44  120.10  63.29 

142.26  104.90  62.77 

138.74  93.02  61.21 

136.57  62.27  60.25 

133.63  58.35  58.96 

132.58  42.74  58.49 


With  EDa  Feeding  air  With  EDa 

100.00  0.0031  0.0431 

62.96  0.0072  0.0812 

47.95  0.0108  0.1095 

41.88  0.0143  0.1390 

37.14  0.0178  0.1674 

24.86  0.0212  0.1875 

23.30  0.0245  0.2078 

17.06  0.0277  0.2232 


3  Li-air  battery  performance  combining  the  effect  of  feeding  air  and  electrolyte 
degradation. 
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discharge  cycle  compared  between  the  cells  with  stabilised  electrolyte  and  with 
electrolyte  degradation  effect.  The  parameters  used  in  the  model  are  the  same  as 
described  in  Fig.  2. 

to  follow  the  straight  line  during  the  cell  simulation  of  8-repeated 
cycle,  then  the  cell  can  estimate  how  many  cycles  will  achieve  until 
the  U2CO3  deposition  show  significant  impact  on  the  Li-air  cycling. 
Fig.  5  also  provides  the  predicted  LhCOa  volume  fraction  extended 
to  the  80th  cycle,  corresponding  to  the  value  of  0.28.  In  other  words, 
when  the  Li-air  battery  is  operated  by  using  ambient  air  with  sta¬ 
bilised  electrolyte  or  without  electrolyte  degradation,  it  can  be 
repeated  discharge  and  charge  to  about  60—65  cycles  before  Li2CC>3 
formation  reaches  the  same  value  as  in  the  case  of  Li-air  battery 
combining  with  electrolyte  degradation  effect  giving  only  8  cycles 
as  compared  in  Fig.  5. 

The  evidence  from  this  study  indicates  that  the  air  feeding  af¬ 
fects  the  cell  performance  in  term  of  decreasing  capacity  due  to  the 
low  diffusion  and  solubility  of  oxygen,  but  this  ambient  air  opera¬ 
tion  does  not  impact  much  for  the  Li2CC>3  generation  during  short¬ 
term  cycles  (8  cycles)  due  to  limited  CO2  contamination  from 
atmosphere. 

- Simulation  at  0.05  mA  cm'2  in  Air 


Specific  capacity  (mAh  gaarbon) 


Fig.  6.  Voltage-capacity  curve  for  a  non-aqueous  Li-air  battery  operated  in  ambient  air 
feeding  on  1st  cycle  validated  with  our  group  experiment  at  a  rate  of  0.05  mA  cnr2. 
The  electrolyte  contains  1  M  LiTFSI  dissolved  in  TEGDME.  The  cell  cycle  is  simulated 
between  2.4  and  4.2  V  versus  Li/Li+  and  both  results  are  tested  at  operating  temper¬ 
ature  298.15  K. 


4.2.  Comparing  the  Li-air  model  with  experiment  using  air  feeding 

As  previous  section  demonstrates  the  Li-air  battery  perfor¬ 
mance  in  various  discharge  rates  on  ambient  air  feeding  without 
being  compared  to  the  real  system,  the  developed  model  in  this 
section  is  now  validated  against  our  group  experimental  data 
owning  to  sufficient  and  accessible  information  on  Li-air  battery. 
The  result  of  1  st  cycle  obtained  from  the  Li-air  battery  assembled  in 
our  lab  is  compared  with  the  simulation  as  shown  in  Fig.  6.  The  air 
electrode  was  prepared  by  loading  3.54  mg  cm-2  of  Super  P 
(Timcal,  surface  area  61  m2  g-1)  without  using  any  catalyst.  The 
electrolyte  loaded  in  the  electrochemical  cell  and  applied  in  the 
simulation  for  this  section  only  was  a  1  M  lithium  bis-tri- 
fluoromethanesulfonyl  imide  (LiTFSI)  in  tetraethylene  glycol 
dimethyl  ether  (TEGDME  or  tetraglyme).  It  can  be  seen  from  the 
graph  that  the  cell  voltage— capacity  curve  at  a  rate  of  0.05  mA  cm-2 
from  the  model  matches  well  to  the  experimental  data  during  cell 
operation  in  1  atm  of  dry  atmospheric  air.  During  discharge,  the  cell 
potential  fell  steeply  at  the  beginning,  from  a  voltage  of  3.1  V  to  a 
plateau  at  around  2.6  V  which  was  in  good  agreement  with  the  Li- 
air  test  cell  for  similar  electrolyte  and  operating  conditions.  The  Li- 
air  model  continuously  discharged  until  reaching  the  cut-off 
voltage  at  2.2  V  corresponding  to  the  discharge  capacity  based  on 
the  weight  of  carbon  alone  at  248  mAh  g^arbon,  whereas  the  real 
system  provided  the  discharge  performance  at  242  mAh  gcarbcm-  The 
validation  of  this  study  indicates  that  the  developed  Li-air  model 
with  the  air  feeding  parameters  can  be  used  as  a  modelling  tool  to 
predict  the  Li-air  battery  performance  in  ambient  condition.  Con¬ 
trary  to  expectations,  the  charging  profile  for  the  Li-air  model  did 
not  fit  well  closely  to  the  experiment  cell  at  the  commencement  of 
the  charging  process  as  demonstrated  in  Fig.  6.  The  former  started 
to  charge  at  4.0  V  and  maintained  at  this  potential  until  the  end  of 
charging  step,  while  the  latter  begin  at  3.6  V  and  gradually 
increased  to  the  same  cut-off  charging  voltage  at  4.3  V.  The  most 
likely  cause  of  different  in  charging  potential  between  the  two  re¬ 
sults  is  that  the  L12O2  formed  on  the  first  discharge  of  the  Li-air 
battery  with  tetraglyme  based  solvent  is  coexisted  with  a  mixture 
of  Li2C03,  HCC^Li,  CH3C02Li,  and  esters,  due  to  electrolyte  decom¬ 
position  [19],  These  products  were  not  included  for  charging  pro¬ 
cess  in  the  case  of  Li-air  battery  model  which  considered  only  the 
Li202  as  the  truly  rechargeable  product.  Flence,  the  charging  po¬ 
tential  of  experimental  cell  may  be  vary  according  to  the  quantity  of 
discharge  products  other  than  Li202.  Moreover,  the  charge  potential 
and  capacity  from  both  experiment  and  simulation  demonstrated 
the  same  amount.  This  comparison  confirms  that  the  model  results 
and  experimental  data  show  similarities  and  the  model  is  a 
promising  tool  to  identify  the  Li-air  battery  mechanisms  and 
forecast  the  cell  performance  in  ambient  air. 

4.3.  Li-air  performance  combining  the  electrolyte  degradation  effect 

In  previous  section,  the  model  demonstrates  only  the  effect  air 
feeding  into  the  Li-air  battery  to  study  the  impact  of  using  air  on 
performance.  However,  in  practical  Li-air  the  challenge  in  non- 
aqueous  electrolyte  degradation  still  remains  during  battery  oper¬ 
ation  and  the  formation  of  Li2C03  is  unavoidable.  Since  the  Li-air 
battery  is  still  in  its  early  developing  period,  there  is  no  electro¬ 
lyte  that  is  perfectly  compatible  to  Li-air  battery  without  electrolyte 
degradation.  Hence,  the  model  is  combined  electrolyte  degradation 
effect  together  with  air  operation.  As  summarised  in  Table  3,  the 
key  results  obtained  from  cycling  Li-air  battery  with  electrolyte 
degradation  are  compared  with  the  cell  using  air  feeding  and  also 
plotted  in  Figs.  4  and  5  in  term  of  capacity  retention  and  Li2C03 
formation,  respectively.  It  is  apparent  that  the  discharge  capacity 
decreases  dramatically  during  cell  cycling  from  250  mAh  gcarbon  in 
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the  1st  cycle  to  very  low  value  of  43  in  the  8th  cycle  due  to  the  high 
amount  of  irreversible  IJ2CO3  deposition  as  present  in  Fig.  5.  This 
cell  performance  falls  faster  than  the  Li-air  with  stable  electrolyte 
as  the  latter  can  maintain  the  discharge  capacity  to  132  mAh  gcarb 
at  the  same  cycle  (see  Table  3).  The  capacity  retention  (decreasing 
to  17%  in  8th  cycle)  from  the  cell  with  electrolyte  decomposition 
also  follows  the  same  trend  of  discharge  capacity  as  shown  in  Fig.  4. 

From  this  study,  the  electrolyte  could  be  considered  as  a  key 
component  and  one  of  the  main  issues  to  be  solved  at  present  to 
sustain  the  rechargeability  of  non-aqueous  Li-air  batteries.  It  needs 
to  stable  in  both  oxygen  rich  electrochemical  condition  during  Li- 
air  operation  and  intermediate-reduced  species  as  well  as  the 
lithium  oxides  (LiOx)  compounds  forming  on  discharge.  Funda¬ 
mentally,  the  stable  electrolyte  determines  a  desired  discharge 
product  (Li202)  and  whether  a  truly  rechargeable  Li-air  battery  can 
be  built  or  not,  i.e.,  high  coulomb  efficiency  during  each  cycle  is 
close  to  99%  to  maintain  battery  cycle  life.  It  is  worth  noting  that 
although  the  electrolyte  degradation  in  this  model  did  not  base  on  a 
detailed  mechanism  which  includes  complex  elementary  steps  in 
the  reaction  sequence  [12,19],  the  generic  Li2C03  formation  reac¬ 
tion  could  be  adequate  to  predict  the  Li-air  deterioration  during 
cycling. 

4.4.  Li-air  performance  with  oxygen-selective  membrane 


operated  in  ambient  air.  As  mentioned  early,  the  immobilized  sili¬ 
cone  oil  membrane  loaded  in  the  various  porous  support  films  have 
been  conducted  in  the  Li-air  batteries  enabling  to  operate  in 
ambient  air  with  better  performance  [17],  To  describe  the  mecha¬ 
nism  of  species  in  membrane,  a  simplified  development  of  the 
theory  of  gas  transport  across  a  membrane  is  presented.  For  the 
single  gas  specie,  several  assumptions  will  be  made  to  simplify  the 
model  and  the  diffusion  of  gas  through  the  membrane  can  be 
defined  by  Fick’s  first  law 


where  N,,m  is  the  diffusion  flux  of  gas  through  the  membrane,  D;im  is 
the  diffusion  coefficient  in  the  membrane  medium,  and  dc,/dx  is  the 
concentration  gradient  of  the  gas  across  the  membrane.  If  the 
diffusion  flux  does  not  change  with  time,  a  steady-state  condition 
exists.  Hence,  the  diffusion  of  gases  through  a  membrane  for  which 
the  concentrations  of  the  diffusing  species  on  both  sides  of  a 
membrane  are  held  constant.  For  a  very  thin  membrane,  if  Dj  m  is 
assumed  to  be  constant,  Eq.  (25)  can  be  rearranged  to 

N,,m  =  Di,m(q  Y'’m)  (26) 


As  one  mechanism  for  U2CO3  formation  and  the  hydrolysis  re¬ 
action  of  the  metallic  lithium  anode  result  from  the  cathode 
exposed  to  C02  and  H20  respectively  from  ambient  air,  it  is 
necessary  to  prevent  these  gases  entering  to  the  Li-air  battery.  For 
example,  the  outer  surface  of  the  porous  cathode  could  be  covered 
with  an  oxygen  diffusion  membrane  that  may  not  only  provide  a 
much  larger  solubility  of  oxygen  concentration  through  the  mem¬ 
brane  but  also  block  C02  and  moisture  content  from  atmosphere 
[16,17],  Moreover,  the  membrane  could  also  minimise  the  evapo¬ 
ration  of  electrolyte  from  the  Li-air  battery.  In  this  section,  the 
developed  model  considers  the  Li-air  performance  with  the  using 
of  membrane  as  both  an  oxygen-selective  medium  to  obtain  a 
better  discharge  capacity  and  a  CO2  barrier  to  decrease  the  Li2C03 
formation.  As  shown  in  Fig.  lb,  the  porous  cathode  of  Li-air  battery 
is  protected  by  the  oxygen  selective  membrane  facing  the  atmo¬ 
sphere.  From  the  previous  section,  the  Li-air  battery  including 
electrolyte  decomposition  effect  tends  to  have  higher  Li2CC>3  for¬ 
mation  than  the  cell  with  the  impact  of  CO2  from  atmosphere  alone. 
Hence,  to  distinguish  these  effects  when  the  Li-air  battery  in¬ 
tegrates  membrane,  the  electrolyte  degradation  effect  is  excluded 
in  this  section. 

A  membrane  could  be  simply  considered  as  an  interphase  bar¬ 
rier  covering  the  Li-air  porous  cathode  and  used  for  the  separation 
of  ambient  air  into  oxygen-enriching  steam  diffusing  through  the 
Li-air  battery.  The  separation  for  each  gas  in  membranes  occurs  due 
to  different  physical  properties  in  term  of  permeability  and  solu¬ 
bility  of  the  species  flowing  through  the  membrane.  In  this  section, 
the  Li-air  model  will  be  integrated  with  the  silicone  oil  as  an 
oxygen-selective  liquid  membrane  for  improving  cell  performance 


The  membrane  characterisation  and  gases  permeability. 

Membrane  properties'1  Value  Ref. 

Thickness  (m)  5  x  KT5  [17] 

Oxygen  permeability  (mol  m~2  s1  Pa1)  2  x  10  7  [43] 

Carbon  dioxide  permeability  (mol  nr2  s_1  Pa-1)  1.08  x  10~6  [43] 

Oxygen  diffusion  coefficient  in  membrane  (m2  s_1)  1.6  x  1(T9  [43] 

Carbon  dioxide  diffusion  coefficient  in  membrane  (m2  s 1)  1.1  x  10  9  [43] 

a  The  membrane  used  in  the  model  consists  of  silicone  oil  as  the  main  material. 


where  c,-  and  cIim  are  the  concentration  of  the  gas  on  upstream  (air 
side)  and  downstream  (after  pass  membrane),  respectively,  and  l  is 
the  thickness  of  the  membrane.  For  the  membrane  characterisa¬ 
tion,  the  air  or  species  permeation  can  be  measured  with  the 
following  equation 


iyim 

Ws 


(27) 


where  F,  is  the  permeation  flow  rate  of  gas  i  (mol  s-1),  Am  is  the 
membrane  surface  area  (m2),  and  Apig  is  the  partial  pressure  dif¬ 
ferential  of  gas  i  (Pa).  By  combining  Eqs.  (26)  and  (27)  together,  we 
can  determine  the  concentration  of  gases  after  passing  the  mem¬ 
brane.  To  date  various  oxygen-selective  membranes  have  been 
developed  and  laminated  on  the  porous  cathode  to  allow  the 
unpressurised  air  permeating  through  the  membranes,  i.e.  no 
pressure  is  applied  at  the  upstream  side  of  the  membrane.  The 
several  assumptions  are  made  to  simplify  the  calculation:  1)  the 
different  partial  pressures  of  individual  gases  between  the  two  side 
of  the  membrane  covered  the  surface  of  Li-air  battery  are  not  know 
unlike  the  membrane-testing  cell  that  can  be  directly  measured  by 
pressure  gauge,  hence  we  assume  this  value  as  low  as  possible,  2) 
the  diffusion  coefficients  of  species  in  the  thin  membrane  are 
assumed  to  be  constant.  The  membrane  characterisation  and  gases 
permeability  for  use  in  the  Li-air  battery  with  oxygen-selective 
membrane  model  are  summarised  in  Table  4. 

The  Li-air  performance  obtained  from  the  air-feeding  cell  inte¬ 
grated  with  the  50  pm  thickness  of  silicone  oil  membrane  is 
compared  with  the  absence  one  as  demonstrated  in  Fig.  7. 
Although,  the  membrane  can  enrich  the  oxygen  concentration  from 
0.62  mol  m-3  in  air  to  2.18  mol  m  3  after  passing  the  membrane, 
this  is  still  not  enough  to  operate  the  Li-air  battery  at  high  discharge 
rate.  Hence,  the  low  discharge  rate  at  0.05  mA  cm-2  is  applied  in 
this  section.  As  shown  in  Fig.  7,  the  silicone  oil  membrane  with  a 
relatively  high  oxygen  permeability  of  2  x  10-7  mol  m-2  s_1  Pa-1 
(oxygen  concentration  2.18  mol  itT3)  at  room  temperature  enabled 
the  Li-air  batteries  with  porous  electrode  to  operate  in  ambient  air 
with  a  specific  capacity  of  796  mAh  gcarbon.  3.5  times  higher  than 
the  case  of  without  using  membrane  due  to  higher  oxygen  con¬ 
centration.  This  discharge  capacity  obtained  from  the  model  is 
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- 0.05  mA  cm"2  in  Air 
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Specific  capacity  (mAh  gaarbon) 

Fig.  7.  Comparison  of  the  voltage-capacity  curve  for  a  non-aqueous  Li-air  battery 
model  operated  in  ambient  air  in  the  case  with  and  without  oxygen-selective  mem¬ 
brane  at  a  rate  of  0.05  mA  cm-2.  The  other  parameters  used  in  the  model  are  the  same 
as  described  in  Fig.  2. 


similar  to  that  obtained  from  the  Li-air  battery  experiment  using 
silicone  oil  liquid  membrane  for  operation  in  ambient  air,  corre¬ 
sponding  to  a  specific  capacity  of  789  mAh  gcarbon  [17].  Moreover,  it 
is  apparent  from  the  graph  that  the  discharge  cell  voltages  also 
increased,  from  ca.  2.60  V  to  2.75  V  with  the  higher  dissolved  ox¬ 
ygen  concentration  in  the  cell  with  membrane  but  the  charging 
voltage  did  not  affect  much  due  to  constant  L^Ch  concentration 
which  is  assumed  as  the  main  product  during  battery  charging.  The 
CO2  concentration  also  increased  after  permeation  through  the 
membrane  owing  to  the  higher  permeability  of  CO2  in  silicone  oil 
than  oxygen  (Table  4).  However,  the  Li2CC>3  generated  from  this 
increasing  CO2  did  not  notably  influence  the  Li-air  performance  on 
the  1st  cycle  when  comparing  the  enhancement  of  specific 
capacity. 


4.5.  Membrane  with  high  oxygen  permeability 


One  of  the  most  significant  properties  of  the  membrane  for 
attainment  in  superior  Li-air  performance  in  term  of  high  capacity 


Oxygen  permeability  x  106  (mol  m'2s 1  Pa'1) 


Fig.  8.  Effect  of  oxygen  permeability  on  maximum  specific  discharge  capacity  at 
different  discharge  rates  of  a  non-aqueous  Li-air  battery  protected  with  a  50  pm  thick 
oxygen-selective  membrane.  The  other  parameters  used  in  the  model  are  the  same  as 
described  in  Fig.  2. 


and  discharge  rate  is  the  high  oxygen  permeability.  In  this  section,  the 
effect  of  oxygen  permeability  in  various  ranges  on  maximum 
discharge  capacity  at  different  rates  is  presented  for  a  50  pm  thick 
membrane.  The  data  were  calculated  using  the  equations  as  derived 
in  previous  section.  This  assumes  that  the  air  window  on  the  mem¬ 
brane  and  the  cathode  electrode  have  equal  geometric  areas.  The 
maximum  specific  discharge  capacity  at  different  rates  for  a  Li-air 
battery  limited  by  oxygen  transport  through  an  oxygen-selective 
membrane  of  any  permeability  can  be  compared  in  Fig.  8.  It  is 
apparent  from  the  graphs  that  the  specific  capacity  enhances  with 
high  oxygen  permeability  and  it  becomes  smaller  at  higher  discharge 
currents  due  to  the  limitation  of  the  oxygen  to  diffuse  inside  the 
cathode.  However,  even  with  the  high  oxygen-permeable  membrane 
the  maximum  discharge  current  to  perform  a  Li-air  battery  is  only 
0.5  mA  cm-2  with  discharge  capacity  below  300  mAh  gcalbon-  The 
results  of  this  study  indicate  that  although  the  specific  capacity  of  Li- 
air  batteries  about  two  times  higher  than  in  conventional  Li-ion 
batteries  (150  mAh  gcalbon).  the  operating  discharge  currents  are 
still  not  high  enough  for  use  in  heavy-duty  power  devices  or  electric 
vehicles  which  require  higher  discharge  rate. 

For  CO2  transport  through  most  oxygen-selective  membranes, 
its  permeability  is  often  comparable  to  oxygen,  i.e.  ratios  of 
Po2 /PC02  are  generally  less  than  unity  [24],  CO2  from  atmosphere 
can  react  with  the  desired  U2O2  discharge  product  or  electrolytes  to 
form  Li2CC>3.  The  former  product  can  reversibly  charge  to  evolve 
oxygen  at  about  4.0— 4.5  V  lower  than  the  latter  which  requires 
high  potential  more  than  4.5  V  to  recharge  [12,42],  Hence,  if  the 
membrane  is  being  used  for  a  rechargeable  Li-air  battery,  it  must 
effectively  separate  the  CO2  from  the  air  inlet  for  maintaining  a 
long-life  cycle  of  Li-air  battery.  This  is  a  challenge  task  because 
materials  with  high  oxygen  permeability  also  have  high  CO2 
permeability  [43,44],  However,  the  CO2  impurity  is  not  the  problem 
in  primary  Li-air  batteries  due  to  the  absence  of  charging  process  so 
the  final  discharge  product  does  not  matter. 

Nitrogen  is  inert  gas  and  has  a  little  or  no  effect  with  lithium- 
based  electrolyte  on  the  Li-air  battery  performance  [24],  There  are 
the  published  studies  describe  Lithium  nitride  (Li3N)  formation 
generated  from  the  direct  reaction  of  nitrogen  with  lithium  at  elec¬ 
trode  in  lithium  ion  batteries  at  room  temperature  during  charge- 
discharge  cycles  [45,46],  The  same  behaviour  could  possibly  occur  on 
the  surface  of  metallic  lithium  in  Li-air  battery.  However,  in  this  work 
we  focus  on  the  porous  cathode  which  contributes  to  the  main  part 
of  Li-air  battery.  Moreover,  the  permeability  ratios  between  oxygen 
and  nitrogen  ( Po2/Pn2  )  are  normally  >1  thus  increasing  the  relative 
oxygen  over  nitrogen  in  the  feeding  steam. 

5.  Conclusions 

A  macro-homogeneous  model  was  developed  and  used  to  analyse 
the  capacity  and  cycling  behaviour  of  the  rechargeable  Li-air  battery 
operated  under  ambient  air  conditions.  The  model  uses  a  set  of  gov¬ 
erning  equations  which  describe  species  transport,  charge  and  reac¬ 
tion  kinetics  within  the  battery,  taking  into  consideration  of  by¬ 
product  formation,  electrolyte  decomposition  and  changes  in  sur¬ 
face  area  and  porosity.  The  model  can  accurately  predict  the  capacity 
feature  and  the  detrimental  effect  of  electrolyte  decomposition  and 
Li2C03  formation  on  the  capacity  retention.  The  model  forecasts  a 
significant  influence  of  using  an  oxygen-selective  membrane,  which 
could  lead  to  up  to  4  times  increment  in  specific  capacity.  The  simu¬ 
lating  results  are  in  a  good  agreement  with  the  experimental  data. 
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Nomenclature 

a  specific  interfacial  area  (m2  m-3) 

Am  membrane  surface  area  (m2) 

d  concentration  of  species  i  (mol  m  3) 

D;  diffusion  coefficient  of  species  i  (m2  s_1) 

Dj  efr  effective  diffusion  coefficient  of  species  i  (m2  s_1) 

E  electrode  potential  of  cathode  at  any  state  (V) 

E°  electrode  potential  of  cathode  at  standard  state  (V) 

/  activity  coefficient  of  LiPF6  salt 

F  Faraday’s  constant  (96,485  C  mor1) 

F  permeation  flow  rate  of  gasses  (mol  s_1) 

H  Henry’s  law  constant 

11  current  density  in  the  electrode  phase  (A  m-2) 

12  current  density  in  the  electrolyte  phase  (A  m  2) 

j  interfacial  transfer  current  density  of  reaction  m  (A  m~2) 

k  reaction  rate  constant 

l  thickness  of  membrane  (m) 

La,  Lc,  L  thickness  of  APL,  separator,  and  porous  cathode 
respectively  (m) 

Mi  symbol  for  the  chemical  formula  or  molecular  weight  of 
species  i  (mol  kg-1) 

n  number  of  electrons  transferred  in  the  electrode  reaction 

Nj  molar  flux  of  species  i  (mol  rrr2  s_1) 
p  surface  effect  factor 

p  partial  pressure  of  gases  (Pa) 

P  Permeability  of  gasses  (mol  m-2  s_1  Pa-1) 

Ti  reaction  rate  term  that  accounts  for  electrochemical  and 

chemical  reactions  (mol  m3  s_1) 

R  universal  gas  constant  (8.3143  J  mol-1  K_1) 

Rfiim  electrical  resistivity  across  L12O2  film  formation  (Q  m2) 

Sj  stoichiometric  coefficient  of  species  i  in  electrode  reaction 

t  time  (s) 

t+  transference  number  of  cation  in  electrolyte 

T  temperature  (K) 

Veen  cell  voltage  (V) 

Zi  valence  of  charge  number  of  species  i 

Greek  letters 

a  transfer  coefficient 

/?  symmetry  factor 

e  porosity  or  void  volume  fraction  of  porous  cathode 

es  volume  fraction  of  solid  phase  in  porous  cathode 

T]  surface  or  activated  overpotential  (V) 

k  conductivity  of  electrolyte  (S  m_1) 

Ke ff  effective  conductivity  of  electrolyte  (S  m_1) 

v  number  of  moles  of  ions  into  which  a  mole  of  electrolyte 

dissociates 

v+  numbers  of  moles  of  cations  produced  by  the  dissociation 
of  a  mole  of  electrolyte 

Pi  density  of  a  solid  phase  of  species  i  (kg  m-3) 

a  conductivity  of  the  electrode  (S  m_1) 

ffetr  effective  conductivity  of  the  electrode  (S  m_1) 

01  electric  potential  in  the  electrode  (V) 

02  electric  potential  in  the  electrolyte  (V) 

A0fiim  voltage  drop  across  Li2C>2  film  formation  (V) 

V  differential  operator 

Subscripts  and  superscripts 
o  initial 

1  electrode  phase 

2  electrolyte  phase 


a  anodic 

c  cathodic 

g  gases 

m  electrode  reaction,  solid  species  or  membrane 
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